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Abstract. We analyze a solar flare observed simultaneously in soft and 
hard X-rays by instruments onboard Yohkoh and the Compton Gamma 
Ray Observatory . Assuming a simple one-dimensional coronal loop that 
is heated by field-aligned currents, we solve the energy balance equa- 
tion to derive the DC-electric field strength necessary to explain the ob- 
served soft X-ray emission by current-dissipation. We use the derived 
DC-electric field to predict the number flux of electrons accelerated by 
thermal runaway and compare this prediction with the number flux of 
nonthermal thick- target electrons implied by impulsive phase hard X-ray 
observations. We find that runaway acceleration can account for the large 
flux (0, 10 36 s~ l * ) of nonthermal electrons provided the loop filling factor is 
£ 1CT 3 such that heating and acceleration occur in filamented structures 
within the loop. 


1, Introduction 

Magnetic field line reconnection models predict DC-electric field components 
parallel to the loop magnetic field (Tsuneta 1995). It has been proposed that 
such DC-electric fields can produce nonthermal electrons in solar flares via ther- 
mal runaway acceleration (Tsuneta 1985, Holman 1985). In particular, a DC- 
electric field will accelerate thermal electrons until a steady state current is 
established. Since the collisional drag on the electrons decreases with increas- 
ing velocity, those electrons with velocities above a critical velocity will undergo 
runaway acceleration producing nonthermal electrons that emit hard X-ray radi- 
ation via thick-target interactions (Holman, Kundu, and Kane 1989). Electrons 
below the runaway threshold continue to heat the loop plasma producing soft 
X-ray emission. This paper investigates the energetics of the DC-electric field 
model by examining whether the DC-electric field strength implied by soft X- 
ray heating is sufficient to account for the number flux of nonthermal electrons 
implied by hard X-ray observations. 

In the following section, we derive the DC-electric field strength by solving 
the energy balance equation in a single loop with uniform cross-sectional area. 
As input to this equation, we use simultaneous soft and hard X-ray observa- 
tions from Yohkoh and the Compton Gamma Ray Observatory (CGRO) of a 
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simple loop flare that occurred on 1992 September 6. From the Yohkoh Soft 
X-ray Telescope (SXT), we infer the geometry of the flaring loop. From the 
Yohkoh Bragg Crystal Spectrometer (BCS), we derive the temperature T and 
emission measure EM of the flare source. From the CGRO Burst and Transient 
Spectrometer Experiment (BATSE), we deduce the number flux of accelerated 
electrons. 


2. Method 


For a one-dimensional loop, the energy balance equation is expressed as 


dU 

dt 


dz dz 


a) 


where U = 3 nkT is the thermal energy per unit volume, Q is the total flare heat- 
ing rate, F c ~ —10 ~ 6 T b ^ 2 dT/dz is the Spitzer conductive heat flux, R ^ an 2 T ~ 1 / 2 
is the optically thin cooling rate, and the velocity gradient term is the enthalpy 
flux of convective motions within the loop. For the radiative cooling coeffi- 
cient, we use a = 2.2 x 10~ 19 which is based on calculations by Raymond, Cox, 
and Smith (1976). We shall assume that the density n and temperature T are 
uniform along the loop length. This assumption is valid about 20-30 s after 
heating onset when thermal conduction will have redistributed the heat energy 
throughout the loop, and hydrodynamic motions will have restored approximate 
pressure balance within the loop (Fisher and Hawley 1990). 

We simplify equation (1) by spatially averaging it with respect to the total 
loop volume such that 


UV = QV C - RV + [bnkTv(0) + F e (0)]A, (2) 

where V is the total volume of the loop and V c is the volume of the current- 
heated region. The enthalpy and conductive fluxes vanish at the loop apex and 
are negligible in the chromosphere ( z = 0) relative to the heating and radiative 
cooling terms. Note that the volume of the current-heated region is assumed to 
be less than the total observed loop volume V = 2 AL, where A is the observed 
loop area and L is the loop half-length. From the ratio of these two volumes, we 
define a filling factor / = V c /V. Following Holman (1985), we assume that the 
heat energy generated within each current channel is distributed into the larger 
loop volume by conductive (or convective) transport processes on a timescale 
less than the heating timescale. 

For a loop that is heated uniformly by current-dissipation, the Joule heating 
rate is given by 


Qcurr = nkTv e (E/E D ) 2 ergs cm 3 s \ (3) 

where u e as 3.2 x 10 2 nT -3 / 2 s _1 is the thermal collision frequency (for classical 
resistivity), E is the electric field strength (assumed uniform along the loop 
length), and Ep = 7 x 10 _8 nT _1 volts cm” 1 is the Dreicer field. The Dreicer 
field is the field strength at which all the electrons in the plasma undergo thermal 
runaway. Substituting Q curr into equation (2), we derive the following analytic 
expression for E 
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E = 


Ep 


(U + R) 

fnkT v t 


w 


Given n and T, we solve equation (4) for the variation of E during the flare 
for different values of the filling factor. Given E , we can subsequently compute 
the rate of runaway electrons assuming that runaway acceleration occurs within 
the same current-heated channels. The runaway rate is given by the formula 


yv ru „ ~ ,35m' e (i?0/J5) 3 * ^ 8 exp[-2 1 / 2 (E D /E) 1 / 2 * - (1/4)(E D /E)]V C sT 1 , (5) 

which includes electrons that axe accelerated out of the thermal distribution 
as well as electrons that are scattered into the runaway regime by collisions 
(Kruskal and Bernstein 1964). The runaway electrons that propagate along the 
loop will heat the plasma by Coulomb collisions, producing thick-target hard X- 
ray emission as they impact the chromosphere. Assuming a power-law electron 
energy spectrum with spectral index <5 above a low-energy cutoff E cy the total 
number flux of hard X-ray-producing electrons is given by 

N thick ~ 3 x 10 33 ai ( 7 - 1) 2 B(7 - 1/2, 1/2 )E C ~~< (6) 

where B(x, y) is the beta function, 7 = <5- 1 is the power-law spectral index, and 
ai is the power- law amplitude, respectively, of the emitted hard X-ray spectrum 
(Lin and Hudson 1976). The low-energy cutoff is determined by the critical 
energy above which thermal electrons exceed the frictional force and undergo 
thermal runaway. The critical energy is given by E^u = m e (Eo/ E)v e 2 / 2, where 
v e is the electron thermal velocity (Holman 1985). 


3. Results 

We apply equation (4) to an M3. 3 flare that was observed by Yohhoh and CGRO 

at 09:00 on 1992 September 6 (Zarro, Mariska, and Dennis 1995). Prom least- 

squares fits of synthetic spectra to the BCS Ca XIX (A 3.177 A) resonance line, 

we obtain the variations of T and EM. Prom preflare SXT images, we infer a 

loop half-length L~3x 10 9 * * * * * * * cm, a loop cross-sectional area A ~ 3 x 10 17 cm 2 , 

and a total observed volume V ^ 2 x 10 27 cm 3 . Prom the emission measure and 

filling factor, we derive the loop density n~ (EMffV) { / 2 . Substituting into 

equation (4), we compute E for different values of /. 

Figure 1 shows the time variation of E for / = 0.1, 0.01, and 0.001. The 
DC-electric field strength increases with decreasing filling factor. This effect 

occurs because the density within the current channels is inversely proportional 

to /. Consequently, since U ~ n, a higher electric field strength is required to 

sustain the enhanced heating rate within the channels. Similarly, the electric 

field strength increases with time during the flare rise phase because the density 
within the current channels increases as plasma is heated and evaporated into 
the loop. 

We use the BATSE LAD Continuous data to derive the nonthormal param- 
eters of the flare. This data type produces hard X-ray spectra in 16 channels 
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Figure 1 . Upper Panel: Variation of DC-electric field strength for 
filling factors / = 0.1 (solid), 0.01 (dash), and 0.001 (dash-dot). Lower 
Panel: Fitted power-law amplitude of nonthermal hard X-ray emission. 


>e tween 10 and > 1000 keV at 2.048 sec temporal resolution. The spectra show 
l well-defined power-law distribution of the form I = in the 20-100 keV 

ange. From a least-squares deconvolution of the LAD response function, we 
>btain the temporal variations of the amplitude a\ and spectral index 7 of the 
>ower-law component. The variation of a* is compared with E in Figure 1 . The 
emporal variations of these quantities are correlated such that the maximum 
)G-electric field coincides with the peak of hard X-ray burst emission. 

For each filling factor, we use the derived values of E to compute the number 
lux of nonthermal electrons from the thick-target equation ( 6 ) and compare it 
vith the number flux predicted by the runaway formula ( 5 ). Figure 2 shows 
Ins comparison for / = 0.001. The latter filling factor yields the optimum 
igreeinent between the measured and predicted fluxes during the rise phase of 
mpulsive hard X-rays. For this filling factor, the peak DC-field strength is 
* * 10 4 volts cm 1 giving a peak runaway flux of 3 x 10 37 s“ l . Such a high 
lux is a direct consequence of the low cutoff energy ( E c « 10 keV) implied by 
he critical energy for thermal runaway. 

For / > 0.001, runaway acceleration alone does not produce sufficient non- 
hennal electrons to match the observed flux. In this case, the density within 
lit- current-heated region is too low to provide a large enough population of 
iiennal electrons to undergo runaway. For / < 0.001, runaway acceleration also 
tils to match the observed number flux. Examination of equation ( 5 ), show* 
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Figure 2. Comparison between the number flux of nontherrnal elec- 
trons implied by hard X-ray observations (solid) and the number flux 
predicted by runaway acceleration for a volume filling factor of 0.001 
(dashed). 

that Nrun ~ exp(-Eo/E). Since Ed ~ n, the runaway rate drops exponentially 
with increasing density. Physically, the runaway electrons become thermalized 
by collisions when the density in the current-heated region becomes very large. 

The predicted and observed electron fluxes show a large discrepancy after 
hard X-ray maximum (as 09:04 UT) and cannot be reconciled with any value of 
/. This discrepancy reflects a breakdown in many of the simplifying assump- 
tions that underly the energy-balance analysis. First, the assumption of a filling 
factor that is constant in time is likely to become invalid as the flare energy 
is distributed throughout the loop system and the heating extends to possibly 
multiple loops. Second, the assumption that the hard X-ray emission is pre- 
dominantly nontherrnal (and thick-target in origin) becomes questionable when 
the temperature within the heated region becomes very high and the apparent 
power-law spectrum is actually the result of a strong thermal bremsstrahlung 
component. The present soft and hard X-ray observations cannot distinguish 
this case. 


4. Conclusions 

Wc conclude from energy balance that DC-electric field heating and accelera 
bon cun self-consistently explain thermal soft X-ray emission and nontherma! 


ri acceleration in solar flares provided that heating and acceleration occur 
•inented subregions with a coronal loop. The fi lamentation is necessary 
u, usure a sufficiently high density of thermal electrons to undergo runaway 
acceleration. Such filamentation of loop plasma is in fact necessitated by elec- 
trodynamic arguments. In particular, for a loop current system to remain stable, 
the induction magnetic field of the current-carrying electrons (including accel- 
erated electrons) within each channel must be less than the ambient magnetic 
field strength. For a cylindrical geometry, this constraint enforces the following 
upper limit on the volume of the acceleration region (Holman 1985) 


Vj £ 


/ c \ 2 / Ep \ 2 
47rn 2 e 2 'v e / V E ' ’ 



where B is the loop magnetic field strength. To satisfy this constraint, the accel- 
eration region must be fragmented into n c = V c /V } multiple filaments. Adopting 
a typical loop field strength of B = 100 Gauss, and using the value of E com- 
puted at hard X-ray maximum for / = 0.001, we deduce a peak n c ~ 10 12 . The 
existence of such a large number of current filaments introduces serious prob- 
lems for maintaining charge-neutrality in the loop plasma. One possible solution 
is that the current systems are closed by cross-field drifting of protons at the 
chromospheric footpoints (Emslie and Henoux 1995). 
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